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Abstract 

The locahzation length of a low energy tightly bound electron-hole pair (ex- 
citons) is calculated by exact diagonalization for small interacting disordered 
systems. The exciton localization length (which corresponds to the thermal 
electronic conductance) is strongly enhanced by electron-electron interactions, 
while the localization length (pertaining to the charge conductance) is only 
slightly enhanced. This shows that the two particle delocalization mecha- 
nism widely discussed for the electron pair case is more efficient close to the 
Fermi energy for an electron-hole pair. The relevance to experiment is also 

discussed. 
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Interacting electrons in random potentials are in the center of recent studies in mesoscopic 
systems because of their relevance to many different phenomena such as the Coulomb gap 
||1|], persistent currents 0, the two-dimensional (2D) MIT (metal- insulator transition) 0, 
and delocalization of two particle states 0,H-§] . 

Although the subject of two particle state delocalization due to electron-electron (e-e) 
interaction has been the center of many recent studies 0, most of those studies concentrate 
on the situation of two particles in an empty band of a one-dimensional (ID) system. This 
is an idealized situation which is not very relevant to experimental systems. A more relevant 
situation, i.e. two excited electrons above the Fermi sea, was considered in Refs. ^j^. In 
Ref. [0 it was shown that two electrons close to the Fermi energy show no enhancement in 
ID systems. Nevertheless, it has been emphasized that the situation may be more favorable 
for higher dimensions 0, for which, up to now, only the empty band case was studied [§. 
Recently it has been pointed out that when the effect of the Coulomb gap is taken into 
account the two electron delocalization scenario becomes less likely. 

In this letter we consider a different type of two particle delocalization, namely the 
delocalization of a tightly bound electron- hole pair (strongly bound excitons). This exciton 
is analogous to the classical compact electron- hole pair excitations discussed in Ref. [|I|, 
which can transfer energy across the system but not charge. This is different from most 
previous discussion of two particle delocalization which were concerned with the charge 
conductance of the system. Thus, delocalization of the electron-hole will manifest itself in 
the enhancement of the heat conductance of the disordered system, an effect which can be 
measured in the laboratory. Indeed, the relevance of the two particle delocalization scenario 
to electron- hole pairs has been noted by Imry 0. In Ref. it was suggested that not only is 
the two particle localization length enhancement relevant to an electron-hole pair, it is also 
much more effective than for the two electron case, because of the Coulomb enhancement 
of the exciton density. This scenario will be verified in this letter by explicit calculation of 
the localization length. Moreover, the existence of delocalized excitons may have interesting 
consequences for low temperature hopping conductance, which will be discussed in more 



details after presenting the evidence for the exciton locahzation length enhancement. 

Our study is based on the following interacting many-particle tight-binding Hamiltonian: 

H = Y. ^kWkjakj - Vj2i(^l,j+i(^k,j + 4+ij«fc,i) + h.c + Hint, (1) 

k,j k,j 

where e^j is the energy of a site {k,j), chosen randomly between —W/2 and W/2 with uni- 
form probability, and V^ is a constant hopping matrix element. The interaction Hamiltonian 
is given by: 

H.nt = U Y: fj^'^'fh'''' (2) 

k,j>l,p \''"k,j ^l,p\/'^ 

where U = e^/b and b is the lattice unit. 

We consider L^ x Ly samples of size 3x2, 3x3, 3x4 and 3x5 with M = 6, 9, 12, 15 sites 
and A^ = 2, 3, 4, 5 electrons correspondingly (i.e., 1/3 filling). The calculation of the electron 
transmission required also the calculation of the A^ + 1 electron Hamiltonian for each case. 
The m X m Hamiltonian (where m = (n)) rnatrix is numerically diagonalized and all the 
eigenvectors l^&f^) and eigenvalues E[^ are obtained. The strength U of the interaction is 
varied between — 30V^ and the disorder strength is chosen to he W = 30V^. The results are 
averaged over 6000, 1000, 200, 25 different realizations of disorder for the M = 6, 9, 12, 15 
samples. 

The zero temperature tunneling density of state (TDOS) for the A^ + 1 electron is defined 
in the usual way pO|] : 

He) = j:\(^r'\T.<^UK)M^ - i^r' - e^)) (3) 

I k,j 

which for the non-interacting case will reproduce the single particle density of states above 
the Fermi energy. The averaged TDOS {i^{e)) = {OAA'y^ f^'^^l-^^! {u{e'))de' (here (. . .) 
stands for an average over realizations of disorder) for different values of interaction is 
presented in Fig. |I]a. Since the distance to the many particle band edge B increases as 
function of the interaction strength U we rescale the single particle level spacing A to offset 
this trend. The rescaled "single particle" level spacing is defined as A' = B/N{M — N), 



which is equal to A in the non-interacting case. For large values A' there is little influence 
of the interaction strength on z/(£:), while for small values of A' there is a change. While 
for U = the TDOS is rather monotonous except for a remanent of the level repulsion at 
A' = 0, there is a signature of the Coulomb gap [|I| at f/ = 301^. Thus, the TDOS shows 
the expected behavior as function of the interaction strength. 

We shall define a tightly bound exciton as the movement of an electron from its ground 
state position to a neighboring site, which may be expressed by the application of a creation 
operator 4j = Jkf + Jkf + ^kf + Jkf (where J^f = alj^^akj and J^f = al^^jUkj) on 
the ground state eigenvector. The overlap of an excited state with a state obtained out of the 
ground state by the creation of an exciton at site (A;, j) is then given by \{'^f^^\Jlj\'^o^^)\'^- 
One may also define the exciton density of state (EDOS) 

Me) = E K^r^l E jU<-'')\'5ie - (Er' - E^'))) (4) 

I k,j 

in an analogous way to the TDOS. The EDOS then describes the weight of a single tightly 
bound excitons at a given energy. The average EDOS {ve{£)) (calculated in the same way 
as the average TDOS) for different values of U is plotted in Fig. |I]b. It is clear that the 
interaction strength influences the EDOS in a different way than the TDOS. There is no 
sign of of the Coulomb gap, instead there is a monotonous enhancement of the EDOS, as 
one may expect from the classical theory of tight bound excitons which predicts a constant 
EDOS as function of e 0. 

The localization properties of a system manifest themselves in the transport behavior of 
the system. For a single particle system one may connect the localization length ^ to the 
transmission ^(r', r ', e) = {r\{e — H + iri)~^\f ') via [^ i{e)^^ = (In |t(r, r ',e)|)/|r'— r '|. 
In Ref. [|12| a many particle Landauer formula connecting transmission and conductance for 
an interacting system coupled to external leeds was developed. For the many particle case 
the transmission in the interacting region is given by 

tifr,^-^ ('■'r'i4ii'.")('''.^i°^-ii'r') ,s, 



where 77 — i> In the numerical calculation we assumed r*and r ' to correspond to the opposite 
corners of the sample, i.e., r = (1,1) and r ' = {l,Ly). 

The heat transport (namely the exciton transport) can be probed for example by attach- 
ing leads which excite a tightly bound exciton at one point of the sample and absorbs it on 
another point. In this case the heat transmission is formulated as 

The localization length of the system are extracted from the numerical data in the follow- 
ing way: We calculate the average log of the (heat) transmission between the two corners of 
the sample for the low lying excitations (In \t{^h){^)\) = (li^((2A')~^| //I^a' ^{/i)(^) ^'j ^O'^^'I))- 
The results are presented in Figs. |^ for (ln|t(ft,)(£: = A')|) and in |pD for ((ln|t(/j)(£ = 
2A')|)). From the slope of (Int) vs. Ly one can deduce the localization length accord- 
ing to ^ oc Ly/ {In t) [|ll|. For the charge transmission we obtain the following localiza- 
tion length ^ = 0.856, 0.866, 0.706, 0.716 at f/ = 0, lOV, 20\/, 30V for e = A' and 
^ = 1.016, 1.206, 1.086, 0.956 at f/ = 0, lOV, 20\/, 30\/ for e = 2A'. Thus, the charge 
transmission is somewhat enhanced for weak interactions, but is suppressed for stronger 
interactions. This is in line with our expectations on the influence of the interaction on 
the charge conductance [0,0. On the other hand for the heat transmission the exciton 



localization length ^h = 0.336, 0.476, 0.596, 0.596 corresponding to f/ = 0, lOV, 20V, 30V 
for e = A' and ^h = 0.376, 0.556, 0.696, 0.726 corresponding to U = 0, lOV, 20V, 30V for 
e = 2A'. Thus, the heat transmission and the exciton localization length are significantly 
enhanced by the interactions. 

Thus, low energy excitons show a significant enhancement of their localization length, 
while electrons in that region show only a moderate enhancement. This can be understood 
as the result of the fact that for electrons interacting via the Coulomb interaction the joint 
density of states of two electrons drastically decreases at small energies due to the Coulomb 
gap in the one-electron density of states. This leads to a suppression of the enhancement of 
the localization length which strongly depends on the two particle density of states. However, 



for excitons, Coulomb effects increase the density of states making these arguments more 
plausible. This can be clearly seen for the stronger interaction strength [U = 201^, 301/) 
for which ^h are strongly enhanced, while ^ is suppressed compared to the non-interacting 
value. Also from a direct calculation of the two electron transmission, the suppression of 
the two electron localization length close to the Fermi energy is evident. 

An enhancement of C,h will of course manifest itself in electronic thermal conduction 
measurements of disordered systems. For a case in which C,h > C,, which has not been 
achieved for the parameters we have chosen in our model, but can not be out ruled for larger 
systems or at a different parameter range, interesting behavior may emerge. For example, it 
is quite possible that for ^ ^ b the exciton localization length may become infinite so that 
the electronic conductivity is exponentially small while the electronic thermal conductivity 
changes as a power of temperature [l^J^. Another interesting effect in such a regime is the 
possibility that at low temperatures the delocalized excitons may assist electron hopping 
much more effectively than phonons 0,^,|I^. This will result in a prefactor of the variable 
range hopping proportional or even equal to a universal value e^/h. Such prefactors were 
observed in the QHE regime (see references and analysis in Refs. |]TB|,|T^) and recently at 
zero magnetic field |T^ . 



In conclusion, we have shown that the two particle delocalization mechanism close to the 
Fermi energy is much more efficient for an electron- hole pair than for two electrons (or holes) . 
This is a result of the enhancement of the exciton density of states near the Fermi energy 
by Coulomb interaction opposed to the suppression of the electronic density of states. This 
results in a large enhancement of the electronic thermal conductance while the conductance 
is only moderately enhanced. 
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FIG. 1. The average transport density of state (a) and average exciton density of state (b) 
as function of the rescaled single electron level spacing of a 3 x 3 system for different values of 
interaction strength. 
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FIG. 2. The average log of the charge (regular symbols) and heat transmission (heavy symbols) 
as function of system length for different values of interaction strength, (a) e = A' (b) e = 2A' 
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